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Abstract 

Background: Culicoides imicola Kieffer and Culicoides bolitinos Meiswinkel (Diptera: Ceratopogonidae) are botli of 
veterinary importance, being vectors of Sclimallenberg, bluetongue and African liorse sickness (AHS) viruses. Witliin 
Soutli Africa, tliese Culicoides species sliow a marked difference in tlieir abundances according to altitude, with C. 
imicola highly abundant in lower altitudes, but being replaced as the dominant species by C. bolitinos in cooler, 
high-altitude regions. 

Methods: The thermal physiology of field collected adults of each species was determined to evaluate whether it 
could account for differences in their distribution and abundance. Critical thermal maxima (CTmax) and minima 
(CTmin), as well as upper and lower lethal temperatures (ULT and IIJ) were assessed after acclimation temperatures 
of 19°C, 24°C and 29°C. Critical thermal limits were determined using an ecologically relevant rate of temperature 
change of 0.06°C.min~\ 

Results: Significant differences in CTmin and CTmax were found between acclimation temperatures for C. imicola 
and C. bolitinos. In C. bolitinos, the LLT of individuals acclimated at 24°C was significantly improved (LLT5o = -6.0rC) 
compared with those acclimated at the other temperatures (LLT50 = -4°C). Acclimation had a weak (difference in 
LLT50 of only rC) but significant effect on the LLT of C imicola. When CT^jn, CT^ax. LLT and ULT were 
superimposed on daily maximum and minimum temperature records from locations where each tested Culicoides 
species is dominant, it was found that temperatures frequently declined below the CT^jn and LLT of C. imicola at 
the location where C. bolitinos was dominant. 

Conclusions: The distribution and abundance of C. imicola is likely directly constrained by their relatively poor 
tolerance of lower temperatures. Results for C. bolitinos suggest that the adult phase is hardy, and it is hypothesised 
that the thermal biology of other life stages could determine their range. 
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Background 

Biting midges in the genus Culicoides (Diptera: Ceratopo- 
gonidae) are incredibly abundant, small flies (ranging from 
1-4 mm). Female biting midges feed on vertebrate blood, 
and in so doing, act as vectors of up to 66 viruses that 
affect animals including humans [1,2], At least three of 
these viruses, namely bluetongue (BTV), epizootic haemor- 
rhagic disease, and African horse sickness (AHSV) viruses. 
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all belonging within the genus Orbivirus (Reoviridae), 
cause diseases of such international significance in live- 
stock that they have been classified as notifiable to the 
World Organisation for Animal Health (OIE), The rapid 
spread of BTV in Europe, outbreaks of AHSV in Spain 
since 1987 [3], and recent outbreaks of a novel orthobu- 
nyavirus, Schmallenberg virus [4,5], has demonstrated the 
devastating effect of these viruses on naive livestock popu- 
lations. Due to the emergence and over-wintering poten- 
tial of these viruses in Europe, whereby novel European 
Culicoides species have been incriminated as the vector 
[6], there is an increased interest in the vector competence 
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and ecology of Culicoides species [7]. Similarly, in sub- 
Saharan Africa where AHS is endemic and widespread, 
long term economic losses have resulted from export em- 
bargoes on horses as well as direct loss of livestock [8]. 
The population dynamics of Culicoides species are mark- 
edly influenced by environmental conditions [9]. On ac- 
count of these factors, there is a need to understand the 
biotic and abiotic drivers of Culicoides distribution and 
abundance. 

Two Culicoides species have been implicated in the 
transmission of BTV and AHSV in southern Africa, 
namely C. imicola and C. bolitinos, Culicoides (Avaritia) 
imicola Kieffer is one of the most widely distributed mem- 
bers of Culicoides in the world, extending over the African 
continent into southern Europe and eastwards to south 
China [10,11]. Based on its wide geographical distribution 
and host preference for larger mammals, as indicated by 
the high numbers collected near livestock, C. imicola is 
considered the principal vector of AHSV, equine encepha- 
losis virus and BTV in South Africa [12,13] and BTV in 
southern Europe [14,15]. Often found in great densities, 
with single trap captures able to reach 10^ individuals 
[16], they may make up 99% of all Culicoides species col- 
lected near livestock in areas where they are abundant 
[10,17]. The larvae develop in moist soils rich in organic 
matter, often in full sun. In years with above-average rain- 
fall, a 200-fold increase in abundance may be observed 
[16] due to increased availability of breeding sites. In con- 
trast, C. bolitinos Meiswinkel occurs in low numbers 
across sub-Saharan Africa, but can become the dominant 
species in cooler, high altitude regions [16,18] such as the 
Free State province of South Africa, and Lesotho. In South 
Africa, Culicoides bolitinos was also shown to be abundant 
in the winter rainfall region of the Western Cape province 
[19,20], and the dominant Culicoides species, in the ab- 
sence of C. imicola, in the sandy dunefields adjoining Port 
Elizabeth in the Eastern Cape Province [21]. The absence 
of C. imicola at Port Elizabeth in light trap collections 
made at Struisbaai and Alexander Bay on the southern 
and western coastline was attributed to the sandiness of 
the soil [22]. Described in 1989 [11], C. bolitinos was iden- 
tified as a potential vector of AHSV as recently as 1998 
[16]. The larval cycle is completed in bovid dung (buffalo, 
domestic cattle and blue wildebeest) [11]. The abundance 
of C bolitinos is thus affected by the availability of suitable 
dung, with the decomposition of dung potentially acting 
to incubate C. bolitinos larvae. 

Despite the role of Culicoides species as potential vec- 
tors and their increased relevance to epidemiological re- 
search [23] little is known about their basic physiology. 
This is most likely due to their small size, nocturnal be- 
haviour as well as ease of injury during handling, which 
has discouraged or hampered experimental manipulation 
[10]. As noted above, C. imicola and C. bolitinos show 



marked differences in abundance, seemingly according to 
regional climate. It is likely, therefore, that thermal physi- 
ology may be a determining factor in their distribution and 
abundance. Most of the physiological processes of insects 
are temperature-dependent [24] and thermal limits are sig- 
nificant for fitness because they determine the ability of an 
organism to remain active and survive during extreme con- 
ditions [25]. Extending from this tenet, thermal limits form 
a fundamental part of mechanistic niche modelling [26] 
that aid in the generation of distribution maps and risk 
management strategies. These distribution models are in- 
creasingly important with a changing climate, in order to 
assist in predicting future outbreaks or range expansions 
or to assess local changes in dominance of species [27]. 
Some work has been done on the role of temperature on 
the geographical distribution of Culicoides species [28], in- 
cluding C. imicola [29-32], mainly through the use of cli- 
mate envelope models. The effect of temperature on BTV 
replication within both C. imicola and C. bolitinos has also 
been investigated [33,34]. However, the tolerance of these 
species to temperature extremes has never been estab- 
lished, and spatial models are coarse [23]. Knowledge of 
the thermal biology of these species may aid in under- 
standing their population dynamics, particularly their dis- 
tribution and abundance, and may even provide insight 
into the epidemiology of the arboviruses they transmit. 

This study determined the thermal limits of adult C. 
imicola and C. bolitinos, Phenotypic plasticity in measured 
thermal tolerance traits was also determined. Results were 
related to the temperature conditions of locations within 
the known geographical distribution of these two livestock 
associated species. 

Methods 

Collections and study animals 

Rearing Culicoides in the laboratory is extremely difficult 
[10], and due to unsuccessful rearing attempts during pilot 
studies, only the thermal limits of adult, field collected in- 
dividuals were investigated in this study. All midges were 
collected using Onderstepoort light traps during March, 
April and October, 2013. Trap catches of both species 
are highly female biased, so all individuals tested using 
the methods outlined below were female. All C. imicola 
were trapped at ARC- Onderstepoort Veterinary Institute, 
Gauteng, South Africa (ARC-OVI; 25° 39' 2.58" S 28°11' 
8.58" E; 1219 m a.s.L). The Onderstepoort area is relatively 
frost-free with only a few days where the minimum 
temperature will fall below 0°C. Culicoides bolitinos were 
trapped on Koeberg farm near the Clarens Valley, Free 
State, South Africa (28° 0' 0" S, 28° 31' 0.12" E, 1631 m a.s. 
1.). Clarens regularly experiences minimum daily tempera- 
tures below freezing, with occasional snow in winter. Both 
species were transported in cool boxes to the Hatfield 
Campus of the University of Pretoria where they were 
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housed as described by Venter et al [35] and fed a 10% 
sucrose solution (in distilled water) from saturated cot- 
ton wool Adult C imicola and C. bolitinos were kept at 
acclimation temperatures of 19°C, 24°C and 29°C with a 
12 L:12D photoperiod for three days before being sub- 
jected to thermal tolerance assays. In most cases, in- 
sects held at constant temperatures for a period of 
three days will exhibit phenotypic changes in thermal 
tolerance traits [36]. 

Temperature data from meteorological stations close to 
ARC-OVI (Pretoria-Arcadia, 30164, -25° 44' 18.8514" S 
28° 12' 26.3874" E, 1400 m a.s.l) and Clarens (Golden 
Gate-Clarens, 30647, -28° 30' 13.716" S 28° 35' 1.68" E, 
1849 m a.s.l) were obtained for the period 1st January 
2011 until 31st October 2013 from the ARC-Institute for 
Soil, Climate and Water (Arcadia, South Africa). These 
data were used to determine the relevance of the thermal 
limits of both species. 

Critical thermal limits 

A dynamic, ramping temperature protocol was used to de- 
termine the critical thermal minimum (CTmin) and max- 
imum (CTmax) for each species after having been exposed 
to the three acclimation temperatures. It has been re- 
ported that critical thermal limits vary with methodo- 
logical context [37], so an ecologically relevant rate of 
change was determined and used. Little of the micro- 
habitat preference of adult C. imicola is known, except 
their heliophilic and exophilic nature [31]. Consequently, 
a known larval habitat and adult catch site was chosen. 
The appropriate rate of temperature change was deter- 
mined using iButtons with a 0.5°C resolution housed in 
waterproof silicone enclosures (SL50-ACC06, Maxim, 
Dallas, TX, USA). The iButtons were buried just beneath 
the soil surface in two waterlogged areas, in full sun, at 
ARC-OVI. Emergence traps indicated this area to be a 
habitat for C. imicola. The mean daily temperature change 
of this microhabitat was found to be 0.06°Cmin"^. 

To determine critical thermal limits, individual midges 
were placed in an aluminium 8-chamber stage, made for 
use under a dissecting microscope (Additional file 1) and 
connected to a programmable circulating refrigerated 
waterbath (CC-K25, Huber Kaltemaschinenbau, Offen- 
burg, Germany). A 1:1 blend of propylene glycol and 
water was pumped through silicone tubing from the 
waterbath to the stage, through which a number of chan- 
nels had been drilled to permit flow of liquid and the ex- 
change of heat. Biting midges were placed in individual 
0.2 mL plastic microcentrifuge tubes (430-11, WhiteSci, 
Bellville, South Africa) and subjected to a heating or cool- 
ing rate of 0.06°C.min"^ after a ten minute equilibration 
period at the start temperature of 24°C. Sample size for 
each treatment was f2 = 14 (two trials, each with n = 7). 
CTmin was defined as the point where muscle contraction 



ceased (i.e. cold stupor). CT^ax was initially defined as the 
onset of muscle spasms [36]. However, no discernible dif- 
ference between erratic movements due to discomfort or 
muscle spasms could be seen, and thus CT^ax was rede- 
fined as the heat stupor point, when muscle contractions 
ceased. In both cases, midges were gently prodded with a 
horse hair attached to a wooden toothpick to ascertain an 
inability to move. A type-T thermocouple connected to an 
eight-channel temperature data logger (TC-08, Pico Tech- 
nology, St. Neots, UK) was inserted into an empty tube to 
record temperature in the remaining chamber of the stage. 

For each species, data were checked for equality of vari- 
ances (Levenes test) and inspected for normality. In all 
cases, the data met these assumptions. The effect of accli- 
mation on both CTmin and CT^ax were assessed by sub- 
jecting data to a one way analysis of variance. Post-hoc 
multiple comparisons were performed using Fishers LSD 
tests. Analyses were implemented using Statistica v. 7 
(StatSoft, Tulsa, Oklahoma, USA). 

Upper and lower lethal temperatures 

A static temperature, plunging protocol was adopted to 
determine lethal temperatures. Ten to 20 individuals 
from each species and acclimation temperature were 
placed in 2 ml centrifuge tubes, with 5 replicates per test 
temperature. The centrifuge tubes were sealed before be- 
ing placed in a resealable plastic bag that was immersed 
in a programmable circulating refrigerated waterbath. 
The fluid (1:1 blend of propylene glycol and water) in 
the waterbath was held at a set test temperature. The 
midges were exposed to the test temperature for 2 hours, 
after which they were placed in a 70 mm Petri dish lined 
with a moistened 2 cm^ filter paper, and given 6 hours 
to recover. Survival was noted after the recovery period. 
Temperatures tested were: -6°C, -3°C and 0°C for lower 
lethal temperatures (LLT), and 36°C, 39°C and 42°C for 
upper lethal temperatures (ULT). Pilot studies indicated 
no temperature-related mortality above (for LLT) or 
below (for ULT) these temperatures and they were omit- 
ted. Complete mortality for C. bolitinos acclimated at 19 
and 29°C was recorded at a test temperature of -6°C, 
and so tests at lower temperatures were performed on 
these flies. Culicoides bolitinos acclimated at 24°C sur- 
vived below -6°C (100% mortality was recorded at -8°C) 
but this data was not included in analysis of LLT. 

Generalised linear models with binomial distribution 
and logit link function were used to test the effect of test 
temperature, acclimation and the interaction between 
test temperature and acclimation on survival. Probit re- 
gression analyses were used to determine upper and 
lower median lethal temperatures (ULT50 and LLT50). 
All lethal temperature data analyses were performed in 
SPSS Statistics 21 (IBM Corporation, Armonk, New 
York, USA). 
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Results 

Critical thermal limits 

Critical thermal limits for C. imicola responded signifi- 
cantly to acclimation (CTmin: ^2,39 = 23.553; p < 0.0001; 
CTn,ax: F2,39 = 22.2; p < 0.0001) with an increase of 2.10°C 
in CTmin and 1.63°C for CT^ax between acclimation tem- 
peratures of 19 and 29°C (Figure 1). Post hoc comparisons 
indicated that the CT^in of C. imicola from each acclima- 
tion temperature was significantly different from all other 
acclimation temperatures. For CT^^, C, imicola accli- 
mated at 19°C and 24°C did not differ significantly from 
each other but both differed significantly from the individ- 
uals acclimated at 29°C. 

For C. bolitinos, variation in critical thermal limits in re- 
sponse to acclimation was less pronounced than in C. imi- 
cola (Figure 1; CT^in: 137°C; CTmax' 1.30°C), but these 
differences were statistically significant (CTmin: F2,38 = 
4.64907; p = 0.0156; CT^^ax: ^2,39 = 11.3; p = 0.0001). Fishers 
LSD tests revealed that the CTmin between 19°C- and 24°C- 
acclimated individuals did not differ significantly, but both 
differed significantly from C. bolitinos acclimated at 29°C. 
For CTmax) the 24°C treatment was lowest and differed sig- 
nificantly from the 19°C- and 29°C-acclimated individuals, 
which did not differ significantly from each other. 



Lethal temperatures 

The interaction of LLT test temperatures and acclimation 
temperature had a significant effect on survival of both 
species (Table 1, Figure 2 A and B). Culicoides bolitinos ac- 
climated to 24°C had the lowest LLT50 (-6.01°C), whereas 
those acclimated at 29°C had the highest LLT50 (-3.98°C). 
Parameter estimates relative to an acclimation temperature 
of 29°C were inspected and indicated that survival of C. 
bolitinos was significantly higher (p< 0.001) when accli- 
mated at 24°C. LLT of those acclimated at 19°C and 29°C 
were not significantly different from each other. There was 
no significant effect of acclimation on the LLT of C. imi- 
cola. The LLT50 between acclimation treatments varied by 
1°C, with the lowest estimated for 19°C (-4.39°C). 

Interactions between acclimation and temperature 
were significant for the ULT of both species (Table 1) 
Upper lethal temperatures were generally similar for 
both species, (mean ULT50: C. imicola = 40.05°C; C. 
bolitinos = 40.39°C). No significant difference relative to 
29°C was noted for 19°C-acclimated C. bolitinos, 
though individuals acclimated to 24°C did differ signifi- 
cantly from 29°C, with 24°C-acclimated C. bolitinos 
having the highest ULT50 (41.59°C). For C. imicola, par- 
ameter estimates relative to 29°C indicated a significant 
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Figure 1 Effect of acclimation on critical thermal maximum for (A) C. imicola and (B) C bolitinos, and critical thermal minimum of (C) C. 
imicola and (D) C bolitinos. In each case, the mean ± 1 S.E. is shown (n = 14). Means labelled with the same lowercase letter are not significantly 
different from each other at the 0.05 level (Fisher's LSD). 
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Table 1 Effects of temperature and acclimation on the 
thermal limits of Culicoides bolitinos and C. imicola 



Species 


Trait 


Effect 


Wald/ 


d.f. 


P 


C bolitinos 


ULT 


Intercept 


276.072 


1 


<0.0001 






Acclimation 


341.425 


2 


<0.0001 






Temperature 


5.826 


1 


0.016 






Acclimation x Temperature 


521.357 


5 


<0.0001 




LLT 


Intercept 


59.813 


1 


<0.0001 






Acclimation 


37.235 


2 


<0.0001 






Temperature 


128.514 


2 


<0.0001 






Acclimation x Temperature 


107.761 


4 


<0.0001 


C. imicola 


ULT 


Intercept 


4.1 


1 


0.043 






Acclimation 


29.793 


2 


<0.0001 






Temperature 


63.892 


2 


<0.0001 






Acclimation x Temperature 


67.629 


4 


<0.0001 




LLT 


Intercept 


167.244 


1 


<0.0001 






Acclimation 


6.978 


2 


0.031 






Temperature 


0.993 


1 


0.319 






Acclimation x Temperature 


425.118 


4 


<0.0001 



difference between the 19°C and 29°C acclimated flies but 
none between 24°C and 29°C. Culicoides imicola showed 
near-total mortality when exposed to 42°C, with only 53% 
of the 29°C treatment surviving, and 0% of the other treat- 
ments. Culicoides bolitinos showed similar mortality at 42° 
C, with 0% survival of 19°C and 29°C The 24°C-accli- 
mated C. bolitinos, however, had a survivorship of 26.13%, 
with complete mortality documented only when subjected 
to 44°C (Figure 2). 

Discussion 

Temperature is a determining factor in insect life history, 
and thermal tolerance plays a key role in determining in- 
sect distribution [24,38,39]. Critical thermal limits are valu- 
able for determining the relative temperature tolerances of 
insects acclimated to different temperatures, or originating 
from different latitudinal and altitudinal ranges [40]. Typic- 
ally critical thermal minima are more variable in response 
to acclimation than critical thermal maxima, both inter- 
and intra-specifically [41,42]. Interspecific differences in 
our study agree with this notion, with the CT^ax of C, imi- 
cola and C. bolitinos being similar (41.87°C SE ± 0.976, and 
41.56°C SE ± 0.828 respectively), whilst the CT^^in of C. 
imicola (1.61°C SE ± 1.175) is nearly a degree higher than 
that of C. bolitinos (0.65°C SE± 1.117). For the Culicoides 
species tested in this study, expected intra-specific vari- 
ation was found only for C. imicola, with acclimation at 
temperatures from 19 to 29°C leading to a change in CT^in 
of 2.10°C, but only 1.63°C for CTmax- In the case of C. boli- 
tinos, while acclimation had a statistically significant effect 
on both CTmin and CT^ax^ the change in CT^in (1.37°C) 



was similar to that of CT^ax (1.30°C). For C. imicola, 
CTmin changed in a stepwise manner, with approximately 
1°C difference between acclimation temperatures. Culi- 
coides bolitinos showed little variation between the two 
lower acclimation temperatures (<0.5°C), but the 29°C-ac- 
climated midges showed a significant increase in CT^in, 
possibly because 24°C still falls within the optimal per- 
formance breadth for C. bolitinos, whereas 29°C might be 
considered an extreme' temperature that reduces per- 
formance (see Deere & Chown [43]). 

Little difference in the upper lethal temperatures was 
found between C. imicola and C. bolitinos, highlighting 
limited variability in ULT. Similar to critical temperatures, 
it has previously been found that variability in lethal limits 
is greater at a family or generic level, with limited variabil- 
ity between species [44] . Upper thermal limits are also less 
pliable than lower limits, among species, populations, and 
over time (i.e. plasticity responses) [26,42] Our study cor- 
roborates this, with little difference in survival of C. imi- 
cola and C. bolitinos at high temperatures.. However, our 
study also points to a dissimilar response of different spe- 
cies to experimental conditions, and possibly environmen- 
tal variation, as suggested by Chown et al. [37]. Culicoides 
imicola had a clear response to acclimation, with prior ex- 
posure to high temperatures improving heat tolerance in a 
stepwise manner. This was not the case with C. bolitinos 
however, as individuals acclimated to 24°C outperformed 
both those acclimated at 19°C and 29°C. The relatively 
poor performance of C. bolitinos when exposed to high 
temperatures after acclimation at 29°C contrasts to the 
idea of the beneficial acclimation hypothesis, and may in- 
dicate a deleterious acclimation response. The deleterious 
acclimation hypothesis predicts substantial negative per- 
formance implications for individuals held at extreme 
temperatures, both relative to those at median tempera- 
tures and those that typically do not experience such ex- 
tremes in the field [43]. Indeed, although not formally 
tested, C. bolitinos acclimated at 29°C exhibited higher 
mortality in holding cages, with no individuals surviving 
past a week in the laboratory, whereas conspecifics accli- 
mated at 19°C and 24°C were capable of surviving in the 
laboratory for at least two weeks. Determining the degree 
to which stress proteins and other thermal protectants are 
up-regulated would provide clarification as to whether a 
temperature of 29°C is indeed detrimental [40], as would 
longevity and other performance assays. 

Acclimation had no significant effect on the LLT of 
C. imicola, with LLT50 values similar for all acclimation 
temperatures (c. -4°C), and total mortality at -6°C. In 
contrast, acclimation had a significant effect on C. boliti- 
nos, with the LLT50 for the 24°C-acclimated individuals 
reaching -6.0 1°C, but no significant difference between 
19°C and 29°C was found. The deleterious acclimation hy- 
pothesis might once again explain the observed pattern, as 
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Figure 2 Effect of test and acclimation temperature on survival of cohorts of C imicola and C bolitinos. Survival was assessed 6 hours after 
each cohort was exposed to the test temperature for 2 hours. (A) Survival at high temperatures by C. imicola; (B) survival at high temperatures by C. 
bolitinos; (C) survival at low temperatures by C. imicola; (D) survival at low temperatures by C. bolitinos. In each case, the mean ± 1 S.E. are shown (n = 5) 
for each test and acclimation temperature combination. 



would a mismatch between optimum temperatures and 
realised habitat. Another possibility may be the uncon- 
trolled age and physiological state of the field-collected in- 
sects. To obtain sufficient numbers for all experiments it 
was necessary to collect additional individuals from the 
field by trapping on several occasions. Senescence has 
been suggested to reduce stress tolerance in insects [41]. 
This is supported by Luckinbill [44], who found longevity 
selection improved low temperature tolerance in Drosoph- 
ila melanogaster [41]. Adult C. bolitinos, however, have an 
average life expectancy of approximately 20 days, but may 



survive up to 90 days [15]. Similarly, C. imicola adults can 
survive for more than 15 days and C. pycnostictus for 
54 days at temperatures of -1.5°C [17]. Given that our col- 
lection methods were consistent and that the period of 
time that elapsed before testing was an equivalent propor- 
tion of the life expectancy of both species, we believe that 
our data are a true representation of the ability of C. imi- 
cola and C bolitinos to withstand temperature extremes. 

This study indicates that lower temperature limits may 
play a role in constraining the distribution of C. imicola. 
Modelling attempts often rely on a mean temperature 
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threshold of 12.5°C to predict presence or absence of C. 
imicola [29,45] as established by Purse et al [46]. Peters 
et al. [45] on the other hand, found absence indicated by a 
mean temperature of 14.4°C and presence by a mean of 
15.7°C. For the two years of weather data we analysed, 
Clarens had a mean annual temperature of 14.5 1°C, two 
degrees above the suggested thermal intercept suggested by 
Purse, but very close to that of Peters and colleagues. Fur- 
thermore, the suggested mean minimum temperature 
thresholds of Peters et al. (Absence: 9.4°C; Presence: 10.0° 
C) coincide strongly with the mean minimum temperatures 
of our study sites (Clarens: 8.67°C; Onderstepoort: 10.74°C). 
However, we propose that mean temperatures are only an 
indication of potential distribution, and it is extreme mini- 
mum temperatures reached that determine population per- 
sistence of C. imicola. The lowest temperature recorded for 
Clarens in 2012 was -8°C (Figure 3), with temperatures 
regularly below 0°C and often dropping below the LLT of 
adult C. imicola in winter. In tandem with a relatively high 
CTmin, it seems likely that the thermal environment in the 
areas surrounding Clarens may prevent adult C. imicola 
from surviving or establishing in any great abundance, con- 
curring with the negative correlation found between C. imi- 
cola abundance and low winter temperatures [47]. In 
addition, it is highly likely that the low winter temperatures 



of Clarens and similar areas affect larval development and 
survival because the lower developmental threshold of lar- 
vae is ~11.75°C (calculated from [10]). Indeed, Baylis et al. 
[30] suggested that minimum temperatures determine the 
ability of the C. imicola to successfully overwinter as larvae 
(for one or two months [47]), as the larvae inhabit the O 
horizon and are easily affected by ambient temperature in 
the open habitats they prefer [31]. Furthermore, Legg et al. 
[48] suggested that median lower developmental thresh- 
olds are useful for determining phenology (and thus also 
distribution potential). 

In the case of C. bolitinos it appears that the thermal 
limits of adults are not related to their distribution be- 
cause their upper temperature tolerance is similar to that 
of C. imicola and they are tolerant of lower temperatures. 
It may be that the thermal limits of other life stages, espe- 
cially the larvae, play an important role in limiting the 
abundance and distribution of C. bolitinos. In warmer en- 
vironments where they are less abundant, the exothermic 
decomposition of dung could lead to dung temperatures 
too high for larvae to survive. Evidence for this scenario is 
currently not available, as when we placed iButtons in fresh 
dung, mean dung temperatures were isothermic with mod- 
erate ambient temperature (c. 14°C), making inference dif- 
ficult. Non-lethal thermal stress effects could also explain 
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Figure 3 Daily temperature range for Onderstepoort and Clarens from the end of winter 201 1, to end of winter 2013 (1 August 201 1 - 
31 July 2013). Horizontal lines represent thermal limits of C. imicola and C. bolitinos: ULT50 and CTmax of C. bolitinos and ULT50 of C. imicolo (red); 
ULT50 of C. imicolo (orange); CTmin of C. imicolo (purple); LLT50 of C. imicolo (magenta); CTmin of C. bolitinos (blue); LLT50 of C. bolitinos (light blue). 
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their limited numbers in warmer regions. Scatophaga ster- 
coraria, another dipteran with coprophagous larvae, dis- 
plays temperature-mediated variation in egg size as a 
maternal physiological response [49] . In this study, Blanck- 
enhorn [49] found reduced egg-size and reduced survivor- 
ship of offspring even at intermediate temperatures of 24°C. 
Similar stress effects could suppress populations of C. 
bolitinos in hotter environments. Modelling done on an 
allopatric, identical sister species C. brevitarsis in Australia 
found that mean maximum temperatures could predict 
BTV infection (and thus C brevitarsis distribution) with an 
accuracy of -75% [50]. Considering the notable similarity 
of C. bolitinos and C brevitarsis in morphology, habitat and 
habit, it is possible that it could be similarly employed. 

Conclusion 

It appears that although factors such as precipitation and 
availability of breeding sites may affect C. imicola distribu- 
tion and abundance, thermal limits play an important role 
in determining abundance and distribution of this species. 
Temperatures lower than the CT^in of C. imicola may pro- 
hibit its dominance in high-altitude areas or its expansion 
into more temperate climes. The presence of C. imicola in 
Clarens, even though at low abundance (<1% of all Culi- 
coides spp. captured) [16], indicates that some individuals 
of this species may utilise protected microclimates to es- 
cape severe cold, similar to some European species that 
overwinter in stables [51]. Another possibility may be that 
C imicola reinvades these areas annually, as Culicoides spe- 
cies have been known to travel long distances by wind- 
mediated dispersal [52]. The results for C bolitinos are 
more complex, as adults of this species show a relative wide 
thermal performance breadth in comparison to C. imicola. 
Further research on the thermal biology of juvenile life 
stages is required to explain the more temperate distribu- 
tion of C. bolitinos. Due to difficulty in rearing Culicoides in 
the laboratory, this resolution may prove to be logistically 
challenging. This study shows that thermal limits and toler- 
ance varies between closely related species and that this 
may play a regulatory role in the presence and abundance 
of species. A better understanding of thermal phenomena 
may contribute to the creation of more realistic risk maps 
and the elucidation of the influence of climate change on 
the presence and expansion of viral diseases. 

Additional file 



Additional file 1: Technical drawing of water-jacketed eight- 
chamber stage, made for use under microscope. Each chamber can 
accommodate a 0.2 ml micro-centrifuge tube. 
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